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ABSTRACT 

We present a mechanism based on internal self-annihilation of dark matter accreted from the galactic 
halo in the inner regions of neutron stars that may trigger the full or partial conversion into a quark 
star. We explain how this effect may induce a gamma ray burst that could be classified as short, 
according to the usual definition based on time duration of the prompt gamma-ray emission. This 
mechanism differs in many aspects from the most discussed scenario associating short gamma-ray 
bursts with compact object binary mergers. We list possible observational signatures that should help 
to distinguish between these two possible classes of progenitors. 
Subject headings: short gamma ray burst, dark matter, neutron star, quark star. 

1. INTRODUCTION 



^ Gamma ray bursts (GRBs) are inter esting highly energetic phenome na for which there is still no definitive explanation 
. on the originating mechanism (see e.g. lPiranl2004HGehrels et al.ll2009l ). Regarding the initial event that triggers a GRB, 
the combination of the observed short variability timescales and huge radiated energies points towards cataclysmic 
Q ' events leading to the formation of a stellar compact object with mass M and radius R releasing a gravitational 
5-H , energy Ai? w GK'P /R « 10"'''^ — 10^^ erg. Specifically, according to their time duration, GRBs can b e classified as long 
C/3 ■ (LGR Bs) if the duration of the gamma-ray signal is larger than 2 s and short (SGRB s ) in other cases (iKouveliotou et ahl 
d \ 119931 ). although the boundary is not sharply delimited (see e.g. iZhang et al.l 120091 : iBromberg et al.ll2012l ). There are 
several pieces of evidenc e in favour of an as sociation of LGRBs with massive stars, especially the facts that they occur 
in star- forming regions (iBloom et all 120021 ) and that supernovae have been foun d in association with a few of them 
(see e.g. iStanek et al.ll2003D . This leads to the collapsar model (|Wooslevl Il993f ) where a LGRB is associated with 
the gravitational collapse of a massive, Wolf-Rayct, star leading to the formation of an accreting stellar-mass black 



. hole (BH). An alternative scenario where the central engine formed after the gr a vitational collapse is a young, rapidly 
\^ ■ rotati ng (ms period) neutron star (NS) has also been considered bv lUsovl ()1992[ ) . [ThompsonI (I1994D andlMetzger et al.l 
\ (|2011h . On the other hand, the identification of the progenitors of SGRBs is more uncertain |NakaJ l2007f ) . This is 
• . partially due to the fact that the detection of SGRB afterglows is technically more difficult, and therefore rarer, than 
' for LGRBs. So far the re is only a handful of detected afterglows, since the first detection in 2005 jGehrel s ct a l. 20051 
S ' IVillasenor et"all 120051) . SGRBs can occur in all types of galaxies, either of a star-forming type or not (jBergeill20l3 ). 
" I , This is consistent with the most popular scenario, that is the merger of a NS-I-NS or NS-I-BH binary system. 
. . • As a main motivation for the present contribution and due to the current uncertainties on the origin of (especially 
J> ] short) GRBs, it is worth investigating other possible classes of cataclysmic events where astrophysical stellar compact 
objects are involved. We will focus on NSs. A NS can be partially characterized by a set of static observables, namely 
mass, ^/ns- radius, i?NS- The structure of a NS has been discussed thoroughly in the past (for a review see e.g. 
[Glcndenn ing 200Q ). It can be briefiy described by a dense core, where matter is in the form of a liquid of nucleons 
(or additional particles inclu ding strangeness), an d a crust where matter undergoes a liquid-solid transition into a 
clusterized myriad of nuclei ([Horowitz et al.ir2005[ ). In more detail, the so-called neutron drip (ND) line signals the 
boundary between the outer and inner crust where neutrons start to leak out from nuclei. The mass in the crust is 
small (up to a few percent) as compared to the whole NS and therefore it could be accelerated, if ejected by large 
energies, to large velocities or even high Lorentz factors. 

The equation of state (EoS) that governs the description of matter in the inside of NSs is currently poorly known. 
Whether matter remains in the form of baryons or deconfined quarks is a matter of debate. The problem of the 
internal constituents in a NS has been thoroughly a nd periodically revisited from the early 1930's when Landau 
predicted the existence of such objects (|Landaul 1X9321 ) from the basis of the existence of neutral hadro nic particles. 
From the discovery of the neutron by Chadwick in 1932 to 1984 when it was hypothesized (|Wittenlll984[ ) that a more 
stable type of matter could be possible (formed by roughly equal quantities of uds quarks, arising from weak decay 
of regular ud matter present in nucleonic matter made of protons and neutrons) a wide variety of possible nuclear 
configurations has been proposed. In this latter quark system case, for example, recent attention has bee n devoted to 
lattice calculations showing the possibility of the existence of the ff— dibaryon, with uuddss quark content (jBeane et al.l 
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l2011f ) and experimental input from the LHC is expected on the formation of a quark-gluon plasma. 

A stable type of extended strange matter woul d allow the existence of stars even more compact than NSs, i.e. hybrid 
or pure quark stars (QSs) (jAlcock et al. IfTOSGl ). The transition from a nucleonic type of matter to a quark phase 
is pr e sently not well unde rstood. From the theoretical point of view, some work has been done (jBerezhiani et aH 
I2003t iHorvath et al.l I1992D on the possibility of quark bubble nucleation in a cold system, where the temperature is 
well approximated by T « 0. As this mechanism is not efficient enough, this process may seem unlikely to occur. 
Another series of works have been developed by Ouyed and collaborators which propose the Quark-Nova model in 
which a massive NS converts explosively to a quark star (|Ouved et al.l 1200^ . The transition is obtained mostly via a 
central density increase due to spin-down or accretion forming a conversion front that propagates toward the surface. 
The outcome is ejection of the NSs outermost layers at relativistic speeds. Alternatively, as we explain in this work, 
additional external agents may be able to produce thermal excitations in these dense degenerate matter systems. 

From indications concerning astrophysical and cosmological backgrounds we know that dark matter (DM) is dis- 
tributed inhomogeneously in the universe and, additionally, could be gravitationally accreted onto individual massive 
astrophysical objects. One of the most popular dark matter candidates is considered to be a Majoran a-type weakly 
intera cting particle (WIMP), although other asymmetric dark matter candidates could exist as well (jKaplan et al.l 
|200{| 1. From the experimental point of view, current direct and indirect det ection experiments try to co nstrain its 
properties. For example, dire ct detection earth-based experiments DAMA (DAM A Collaboration! [2010f ). CoGeNT 
(|CoGeNT Collaboratioi3l2011f ). CRESST and lately COUPP (jBehnke et a ll 1201 2, ') seem to point towards a light par- 
ticle in the ~ 10 — 20 GeV/c^ mass range (although this is in tension with other null searches) while there has been a 
lot of interest in indirect detection of new excesses of extra photon components from the G alactic Center that could 
be due to annihilation of a ~ 130 GcV/c^ particle (jBringmann et al"]l2012l : lWenigerll2012| ). albeit so far at only the 
2 — 3(7 level of significance. 

DM passing through astrophysical massive objects is expected to scatter off their constituent nuclei, lose energy, and 
become trapped by their gravitational field. Multiple scatterings cause DM to sink towards the center and afterwards, 
assuming that it may self-annihilate, a series of end products are available. Based on the previous arguments, a set 
of current experiments try to find evidence for this process in the sun, eart h (jAvrorin et al.l[2009t iBaer et al.ll2004[ ) or 
even more compact-sized objects such as NSs (jKouvaris fc TinvakovllmnI) . 

Recently, new ideas on the poss ibility that quark bubbles may be formed by spark seeding in the inner cores of 
dense NSs have been proposed bv iPerez- Garcia et al.l ()2010f ). Sparks would be generated in a sort of Trojan horse 
mechanism by self-annihilation of dark matter particles gravitationally accreted by NSs gathering in their internal 
cores. As estimated, this would allow rel ease energies of up to about one order of magnitude higher than the quark 
binding energies in the present baryons (jPerez-Garcial[2012h . This mechanism then could lead to the conversion of 
NSs into QSs if the transition to quark matter phase is macroscopic. In this way the matter component in the form 
of DM would not only influence cosmological observables, but even NS dynam ics ("Pc rez-Garcia fc Silkll20Tll ). Similar 
arguments have been claimed in the solar seismology observables (jLopes fc Silk,2010l) . 

We investigate in the present paper whether this type of non-repeating cataclysmic event could produce SGRBs and 
what could be a possible specific signature. After a short summary of the scenario in ^ in ^we compare the energy 
released by SGRBs with the energetics of the progenitor model that we propose and in fj4l'^6 discuss the short GRB 
event rate in light of this scenario as compared to observed rates, as well as the natural delay time between the regular 
NS phase and QS formation and how this mechanism may affect the properties of the host galaxies of short GRBs. 
In fJS]we present and analyze the ejecta and crust masses that could be expelled due to the NS— ^-QS conversion and 
we present the details of the resulting central engine model. In ^we attempt to discuss the expected temporal ( ij6.ip 
and spectral ( §6.2p properties of the prompt gamma-ray emission and the afterglow. In ^ 36. 3! we discuss our results 
and possible genuine additional signals and gravitational wave emission that could finally lead to a clear identification 
as compared with the short GRB scenario involving binary mergers. Finally in ^ we give a brief summary of our 
conclusions and suggestions on how to proceed in further work. 

2. SUMMARY OF THE SCENARIO 

There h as been so me previous work where DM accretion in astrophysical objects has been considered () Goldman fc Nussinovl 
[l98 9: Ko uvaris fc T invakov 2010: Lionel et al.,,2012i ). As a summary, we start in this section by motivating the scenario 
where the emission of a SGRB may occur. 

In our previous work (jPerez-Garcia et all 120101 ) the possibility of spark seeding in NS as a result of DM accretion 
was proposed. We found a relationship between the mass of the DM particle, m^^, and the binding energy of the lump 
of quark matter or strangelet formed, E^^^{iii{nA) , mi , B) as calculated from the MIT bag model (jChodos et al.|[T97^ . 

2fm^c^>E±,{MnA).m,,B), (1) 

where ha is the baryonic number density and /i^ and are the chemical potential and mass of the quark of ith— type 
in the light sector, i = m, d, s. B is the bag constant. The channels available to the x annihilation are, so far, not 
known and this uncertainty is phenomenologically parametrized into the fraction /. Inside the NS, it is believed that 
DM scatters one or more times and thermalizes according to a Boltzmann distribution. The thermal radius is given 
by 

/ 3M{t) y/^ 
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Fig. 1. — Qualitative scheme of transition from an initial configuration (A/i,_Ri) to a final one with (A/p,i?p). In blue (red) we show 
difli'erent compact star behaviour for NS (hybrid QS) adapted from[Boinbaci ct al. (2004). Black line depicts the possible transition from 
one configuration to another, which could be triggered by internal self-amiiliilation of DM. 

This radius is, in general, time-dependent since evolution proceeds for central temperatures and baryo nic mass 
densities p^. DM at high density may self-annihilate to channels involving photons, leptons and light qq pairs (jKuhlenI 
I2010D . Due to the fact that either a light or heavy self-annihilating particle in the ^ 1 — lOOGcV/c^ mass-range may 
release a spark of energy able to liberate the constituent quark content of the baryons forming the very dense matter in 
the core of these objects, a pure nucleonic NS could then be considered a mctastable state of compact stars. Note that 
a self-annihilating ~ 20 GcV/c^ particle as recently suggested in the context of direct detection experiments would, in 
particular, fit into this scenario. The annihilation power is then given as 



where is the DM mass density and < a^^v >~ 3 x 10 ^^cm'^/s is the thermally averaged annihilation cross section. 

If the spark seeding is able to form quark lumps or strangelets, most likely off-center, inside a volume ^ r^j^, they could 
individually grow o r coalesce from a stable minimum mass number, so far poorly known, Amin ~ 10— 100, according to 
previous estimates (|Madsenll2001[ ). In this way a burning front may proceed if particular thermodynamical conditions 
are fulfilled. If the burning is partial, the final configuration object would be hybrid, constitued by an inner quark 
matter phase and an outer nucleonic phase. However, if the burning is complete the resultant object would be formed 
by a pure quark phase. In both cases the final configuration in th e mass-radius phase space (MF,i?F), as computed 
from different EoS using a variety of many-body effective theories (jHaensel et al.l[2C)07| ). is a smaller, slightly lighter, 
more compact object than the initial one characterized by {M{,R\). 

In Figure [1] we show a qualitative scheme, inspired from the work by iBombaci et afl (j2004[ ). of how the transition 
would produce a more compact object. The transition point, and therefore the initial configuration mass, would be 
influenced by the presence of spark seeding by self-annihilating DM in this model. 

In such a NS QS transition, part of the outer crust of the initial NS can be expelled, and possibly accelerated 
to relativistic speed, which may lead to a transient episode of high-energy emission. In the following, we discuss a 
possible association of some GRBs with such events. We focus on SGRBs for several reasons: 

• the energetics of bright long GRBs may be a challenge for this scenario, whereas typical energies of short GRBs 
are in a much better agreement (see f^l); 

• a delay of at least A ~ 10'^ yr is expected between the formation of a NS and the transition to a QS, which 
is too long for long GRBs as supernovae within a few days have been found in association with some of these 
events. On the other hand, the scenario naturally leads to a broad distribution of the delay At, allowing long 
delays that are necessary to explain the properties of the afterglows and host galaxies of short GRBs (see 

• the scenario naturally leads to short (or very short) timescales for the gamma-ray emission, pointing towards 
short or ultra-short GRBs (see §6.1|) . 



Since 2005 and the first accurate localization of a SGRB (jGehrels et al. II2005D , the redshift of an increasing number 
of short GRBs has been measured, confirming their cosmological ori gin and leadiiig to a large d ispersion in the isotropic 
equivalent energy released in gamma-rays, E'^^so — 10"^* — 10^^ erg ()Nalaill2007l : lBergeiil2007D . This energy is emitted 
by a relativistic outflow which is probably beamed within an opening angle 6^ so that the true energy release is 




(3) 



3. ENERGETICS 
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= ^E^_iso, where the beaming factor is defined as 

/b=(^) =(l-cos0j)-^^| if0j«l. (4) 

Unfortunately, the op ening an^le 0i cannot be easily measured in GRBs. It can be estimated if a jet break is identified 
in the late afterglow ()Rhoad s"1997^. It is only i n very rare cases that a precise constraint on 9j is obtained, such as 
6'j ~ 7° in GRB 051221 A (ISoderberg et al.l[200l . and even in su ch cases, it remains partially model-dependent. The 
distribution of /b for SGRBs is therefore unknown (INakar|[2007f ). Highly beamed SGRBs with 9j ~ 1° would lead to 
/b = 6.5 X 10'^ whereas larger opening angles of about 9j ~ 30° give /b — 7.5. This will remain as a major source of 
uncertainty in what follows. 

In SGRBs, the isotropic equivalent ene rgy E^,iso is r eleased in gamma-rays on a very short time scale 2 s) 
and observed well above 1 MeV (see e.g. iGuiriec et al.|[2010lV Then, the gamma-ray emission must be produced 
from relativistic ejecta to avoid the compactness problem: this provides a con strain on the minimum value of the 
Lorentz factor for the emission region to be optically thin for 77 an nihilation (|Reej 19661: IBaring &: Harding 119971 : 



ILithwick fc Sarill200lD . The an alysis of a few pre-Fermi short bursts (jGhirlanda et al 



I200I leads to r > Tn,in ^ 15 



from this 77 opacity constraint (|Nakaiil2007[ ). When a SGRB is de tected up to the GeV range by Fermi-L AT , stricter 
constraints are obtained, such as T > 300 - 1000 in GRB 090510 (jAckermann et al.|[2010l: iHascoet et al.ll2012[ ). 

Let us compare these constraints on the energetics and the relativistic nature of the emitting material in SGRBs 
with the predictions of th e scenario described in the previous section. If there is a transition triggered from the inside 
([Perez- Garcia et aI]|2010D . namely due to DM self-annihilation in the central core of a NS (A%g, i?Ns) then the final 
configuration may lead to a more compact object with partially deconfined quark content. The energy gap to this new 
configuration (A/qs, Rqs) will release approximately (here we assume A/ns ~ A/qs for the sake of comparison). 



V ^QS ^NS / 

Note that the intern al energy due to the matter equation of state has been estimated to be of the order of iJgrav 
(|Bombaci et al.|[2000D . therefore this estimate is correct up to a numerical factor that somewhat depends on the EoS. 
The energy Ai^grav liberated in a NS— ?>QS transition can be compared with the energy AE'grav.prog-s-NS — GAf^s/^NS 
released in the gravitational collapse leading to the formation of a NS (standard core-collapse supernova). We define 
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this ratio (in the approximation A/ns ~ Mqs) as /qs/ns ~ I A£;grav,NS^QS / A-/l/grav,prog— >NS I 
calculatio n is made assuming a transition of a hadronic type star with a combined EoS model including additional 
hyperons (jBombaci et al.ll2000l ) to a hybrid star with a quark conten t described by the M IT bag model. Two relativistic 
non-linear Walecka-type EoS models of Glendenning-Moszkowsky (jGlendenningj [19911 ) are tested, namely, GMl (blue 
curve) and GM3 (red curve) usually used to describe the hadronic phase. The typical compressibility values K for 
the GMl model arc stiffer than for the GM3 model. Note that in this scenario of hybrid stars the maximum mass for 
the hybrid object must be somewhat lower than in the case of the pure hadronic star. Although both of them predict 
maximum hadronic stars masses somewhat on the low side with Aff„ax ^ 1.9Af0, (as compared to recent measurements 
reporting masses of about 2A^o ([Demorest et al.ll2010( )). they can be considered illustrative to show a general trend of 
the typical original configurations before transitioning to hybrid or pure QSs. In this figure the value of some nuclear 
parameters such as the s-quark mass has been taken to be nig = 150MeV/c^ and the strong coupling constant as = 0. 
Surface tension is ct = 30MeVfm~^. No curvature energy has been considered. 

In this situation the initial configuration masses and radii are not the same, since they depend on the value of the B 
bag constant in the MIT model. Following iBombaci et al.l ([20000 we consider that the transition takes place between a 

configuration in the initial hadronic star with gravitational mass Mf = PBir)d^r and with some baryonic mass 

M,p = /q'^"^ nB{r)(Pv and that of a final hybrid object (computed with the hybrid EoS) with the same baryonic mass, 

Mp = Afj^, since that charge is conserved. /3_B(r) and nB{r) are the energy density and baryonic mass density for 
the baryonic object. Physically, the transition occurs due to the fact that as the density or central pressure exceeds a 
threshold value, it is energetically allowed that, at the transition point, bubbles of quark matter form. This is due to 
the lowering of the chemical potential in the system as more degrees of freedom are available. These lumps of quark 
matter can grow very rapidly, driving a burning front and the original hadronic star will be converted into a hybrid star 
or quark star. In the mechanism involving self-annihilating DM this energy released could be considered as sparks that 
would coalescence or ignite the medium and allow quark deconfinement out of the baryons present (this is B-dependent 
and therefore, EoS dependent in our model). If we further include corrections due to slightly different masses for the 



initial and final configurations we have /qs/ns ~ |Ai;g,.av,NS^Qs/£'grav,prog^Ns| 



V A/ns ) \ Rqs 



We typically 



find /qs/ns ~ 0-1 ~ 0-3, which illustrates that the transition mechanism discussed here indeed belongs to the class 
of the most energetic phenomena in the Universe. The efficiency /qs/ns is plotted in Fig [5] as a function of the bag 
constant B^/'^. It shows that the efficiency /qs/ns is higher as the bag constant B grows, pointing to heavier DM 
particles in order to ignite the nucleation of quark matter lumps but well in the conservative range shown in Fig. 2 
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Fig. 2. — Efficiency of the energy extraction by the NS— >QS conversion as a function of the bag constant in the combined model of GMl 
(blue curve) and GM3 (red curve) and the MIT bag model. See text for details. 



in iPerez-Garcia et al.l (|201Cl[ ). Further work is needed to clarify the role of the EoS in the efficiency of the transition 
since the possible existence of a burning front able to fully convert the NS into a QS and proceed towards higher radii 
is a crucial point. 

Numerically, using typical values Mns — 1.5M0, _Rns = 12km and i?Qs ~ 7km results in ASgiav — 3.5 x 10^"^ erg. 
In the NS— >QS tr ansition mo st of th e energy Aj^^ra v will be radiated as neutrinos and photons, as showed by detailed 
calculations dJaikumar et al., 2002: Vogt et al.ll2004| ). Additionally gravitational radiation should also be emitted (see 
^6.3|) . We assume now that a small fraction, /ej, of this energy may be injected into the outer crust which is then ejected 
and becomes relativistic. The kinetic energy of the expelled crust is i^kin — /ej AE^gi-av if acceleration is complete. The 
isotropic equivalent energy that could be radiated as gamma-rays by such an ejecta can be estimated by 



^^.a.o - 3.5 X 10- {^^{^^{^^{?:^\ ( -r^ ) erg, (6) 



M 
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where is the efficiency of gamma-ray energy extraction from the ejecta and could range from ~ 0.01 — 0.1 for the 
extra ction of kinetic energy by internal shocks (iRees fc Mesz aros"1994': i Kobavashi et al.lll997t iDaigne fc MochkovitchI 
19981) to ~ 0.5 for p hotospheric emission (see e.g IRees &: Mc szaros 20051) or magnetic reconnection |Th ompson 199^ 
Zhang fc Yan|[20Tl[) . As discussed below, relativistic motion of the ejecta is favored by a small ejected mass, Moj 



10 M0, which corresponds typically to the outer crust. The crust is defined as the region where the density drops 
below nuclear saturation density (NSD) at ^ 2 x 10^"* gcm"'^, and the outer crust is limited by the neutron drip (ND) 
at ~ 4 X lO^^gcm"^ (see Then, the value /oj ~ 10~^ used in Eq. [6] would correspond for instance to a situation 
where 1% of AiSgiav is injected in the crust and 10% of this energy is injected preferentially in the outer crust, so that 
/oj ^ 0.01 X 0.1. As there is la rge uncertainty in this estimate we use a safe approximation below the /cj ~ 0.1 quoted 
in other detailed calculations (jOuved et al.ll2005l ). 

This estimate of i5-y.iso is in reasonable agreement with observations of SGRBs: the NS— S'QS conversion scenario 
investigated here can reproduce measured energies in SGRBs for fhf-yfcj — 3 x 10"^ — 0.3. 

Regarding the relativistic nature of the ejecta, the maximum Lorentz factor, Fmax, that can be reached depends 
on the expelled fraction or ejected mass, Mej, of the crust. The maximum Lorentz factor can be deduced from the 
estimate of i^kin above, and from the ejected mass A/oj as: 

____ E^in ^ [ M \[ Rns]-' [ M Yf y 
" Afejc2 - 10-3 A 12 km j Mq J [ lO'^ Mq J ' 

Again, it seems that Lorentz factors above 15, in agreement with the observational constraints described above, can be 
reached as long as the ejected mass remains low (Mej < lO"** Mq) and the fraction of the energy injected in the outer 
crust is not too small (/ej > 10-^). Even ultra-high relativistic ejecta with F > 100 could in principle be produced, if 
the ejected mass is really small. This could be the scenario where the outer crust is expelled. 

We conclude that the DM self-annihilation triggered NS— >QS scenario can, in principle, release enough energy to 
power a SGRB, and that the ejected crust can reach high Lorentz factors, which is necessary to emit gamma-rays. 
This is however strongly dependent on the two parameters /ej and Mej which are very difficult to predict accurately 
at this stage. This will be briefly discussed in fj5] 
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7^sGRB/ns "4s,max ^ (8 X IQ-* ^ 3 X 10"^) ( ^ j . (10) 
From these estimates only, one can conclude that NS — > QS transition can be much more frequent than SGRBs, 



4. EVENT RATE AND DELAY BETWEEN THE NS FORMATION AND THE TRANSITION TO A QS 

Assuming that the transition of a NS to a QS triggered by the self annihilation of accreted DM in the core can inject 
enough energy into relativistic ejecta to produce a SGRB, it is worth comparing the observed rate of these phenomena 
with the predicted rate of NS— )-QS conversions. The local rate per unit volume of SGRBs, T^-sgrb, can be estimated 
from their observed rate and distribution o f flux using a p opula tion model (luminosity function, comoving rate, etc), 
as has been done bv lGuetta fc PiranI (|2006[ ). iNakar et alj (|2006( ). who obtain 

TesGRB ^ (400 -> 1500) (^^^ Gpc-3 yr-i . (8) 

The lower limit corresponds to a comoving SGRB rate that follows the cosmic star formation rate with a long delay, 
and the upper limit to a constant comoving rate. Unfortunately, the unknown distribution of the beaming factor and 
its average (/b) is again a major source of uncertainty in this estimate. 

In the scenario presented is this work, the local rate of NS formation, estimate d as the local rate of type II supernovae, 
gives an upper limit for the rate of NS — ^ QS conversions. It is of the order of (jPahlen et al.ll2004[) 

<'s''"is.max ^ 5 X 10^ Gpc'^yr"! . (9) 

Then the ratio of the former two rates is 

,(SNii) ~ (8 X 10-4 ^ 3 X 10-3) (^M 

S — > QS transition can 1 

depending on the fraction of NS that will experience such a transition. If all SGRBs are due to NS — >■ QS conversions, 
the low ratio obtained in Eq. [TU] could be (i) either due to the fact that only a fraction of NS — >■ QS transitions lead 
to short GRBs; (ii) or be related to the fact that only a small fraction of NS are converted into QS. 

Having only a small fraction of NS converting to a QS is expected if the delay At between the formation of a NS 
and its conversion to a QS is usually long, e.g. At ~ several Gyr. For a broad probability distribution p{At) of the 
delay At, with a high mean value (At) of the order of the Hubble time, most conversions would not have occured 
yet, resulting in a low ratio 72.sGRB/72-NS^QS,max, the observed short GRBs being produced by the conversions with 
the shortest delays. In addition, such a distribution of p{At) would make the scenario in good agreement with the 
properties of short GRB afterglows an d host galaxies: co ntrary to the case of LGRBs which are always observed in 
central regions of star forming gala xies (iBloom et al.ll2Q0^ , SGRBs can occur at any place (sometimes at the periphery 
or outside) in any type of galaxy (jBergedl2011h . This indicates a delay between the end of the life of the massive 
progenitor star and the production of the short GRB. It is already a well known fac t that most NSs are born with a 
natal velocity kick v that can be as large at « 10^ km/s (jArzoumanian et al.l 12002} ). Therefore, with high values of 
At, not only any correlation with star formation is lost, but the NS can travel a distance 

D^wAi^lOkpcf -)(—^!—] (11) 

^ VlOOOkms-i/ VlOMyry ' ^ ' 

and most likely experience the conversion to a QS far from the galactic central regions, or even outside the galaxy. 

In the scenario presented here a lower limit on At can be obtained from the physics of the transition, which is 
triggered by self-annihila tion of accreted DM from the galactic halo. In our Galaxy the DM density distribution can 
be taken to be of a type (|Navarro et al.|[2"004f ) 

p^(r)=p_2e^[(^'°-^l (12) 

with parameters p_2 ~ 0.22 GeV/cm'^, a — 0.19, and r_2 = 16kpc so that at the solar neighbourhood the Keplerian 
velocity is w ~ 220 k m/s and the local DM density is p^_o ^ 0.3GeV/ cm^. DM from this halo can be accreted by 
gravitational capture (jGouldl[l987D at a rate at peak of NS distribution (jLorimerll2004[ ) of 

C « ^^4^ — particles s'^ , (13) 
TOx(GeV) Px,o 

but off-peak, at about solar circle it may be reduced somewhat. In this estimate, we assume a WIMP-nucleon (spin 
independent) cross section a = u^n = 7 x 10"^^ cm^ (jBertonell2010D . Once the steady state, resulting from competing 
processes of annihilation and accretion, has been reached, the elapsed time is tdm ~ , ^ where V is the 

volume of the star. Typically, if we assume a light ~ 20 GcV/c^ DM particle, as direct detections experiments seem to 
preliminary suggest, we obtain tdm ~ 3.5 x 10"^ yr. For a velocity (temperature) dependent cross section and heavier 
DM particles ^ 100 GeV/c^, this delay can be longer, tdm > 6 x 10^ yr (jKouvarisI I2008D . From the condition that 
the transition should likely occur after reaching the steady state, we get a lower limit on the delay At, 

AW ^ TDM ^ lO-'' - 10^ yr . (14) 
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Notice that this delay is already much too large to allow to consider producing LGRBs wit h NS — > QS conver sions, 
as an associated SN is sometimes found in association with such events within a few days (jStanek et alj 12003'). On 
the other hand, it is too small to have SGRBs uncorrel ated with sta r formation, which is required by the observed 
diversity of SGRB host galaxies and afterglow locations (jBergedl2011h . 

We conclude that At is related to the microphysical processes happening in the burning of the NS. This, in turn, 
depends on the ability of the compact object (mainly related to the EoS) and environment conditions to accrete DM 
and, therefore, cannot be simply predicted and should present some spread and extend up to a large value Atmax that 
can be constrained by observations. 

In our Galaxy, there are hundreds of confirmed NSs, with measured radii and masses (jLattimer &: Prakashll2005() . 
On the other hand, the subsample of confirmed NSs with an estim ate of the age is very small, due to the difficulty of 
the age determinat ion, based on cooling theo ry (jPage et al.ll2009f ). It seems that NSs with an age of at least ^1 — 10 
Myr are identified (jYakovlev fc Pethickll200^ which would put a lower limit on the maximum value of the delay At, 

At,„ax > lOMyr. (15) 

Such delays, as shown by Eq. [11] lead to travelled distance of the order of 10 kpc or more, reaching the outskirts 
of a galaxy but are still too small to de-correlate from the star formation activity. One can not of course exclude 
that much older NSs are present in the sample of confirmed NSs, which would increase At^ax up to 100 My or more. 
Nevertheless, a detailed comparison with the properties of SGRB host would require the knowledge of the distribution 
p(At), which seems yet out of reach. 

Another possible approach to constrain the rate of NS— >QS transitions in the scenario proposed here is to focus 
on the kinematics of observed pulsars. The pulsar tangential v elocity data can be fitted with a bimodal distribution 
peaked around vi = 300 km/s and an upper V2 = 700 km/s (jArzoumanian et al.ll20d^ . The high velocity pulsars 
with V > V2 roughly account for the 10% of the pulsar population. The first peak in the distribution is believed to 
be due to the kick ve l ocity given to the NS when it is formed in a core collapse supernova. It has been suggested 
bv iBombaci &: Popovl (|2004[ ) that the sec ond higher velocity peak could be due to a second kick when the NS— >-QS 
conversion takes place. In previous work (jPerez- Garcia et al.l[2?)T(il : IPerez-Garcia fc Silkll20TTh . it was shown that DM 
seeding in NSs may form a stable and long-lived strange quark matter (SQM) lump that could induce a partial NS 
conversion into a hybrid SQM star. This event may produce, as a consequence, an effect on the kinematics with a birth 
kick and rotation as a result of the partial burning of the NS (jPerez-Garcia fc Silkll201lh able to produce up to u w 10^ 
km/s and relative changes in the angular velocity of Ar2/i7 « 10~^ — 10~^ as a result of the off-center mechanisn i 
driving the transition. Other mechanisms quoted in the literature rely on some kind of asymetry (jLugones et al.ir20Q2D 
in the front progress. From this interpretation of the observed pulsar velocity distribution, ^ 10% could be taken as 
an upper limit for the frequency of the NS— J-QS transitions, leading to an upper limit for the NS— >QS rate which is 
more constraining than in Eq. IIOI i.e. 

^^tiQS,max^5xl04Gpc-V-^ (16) 

Since this rate remains higher than the SGRB rate, at this stage the only possibility in this scenario is to to assume 
that only a fraction /sgrb ~ 1 — 10% of NS — >■ QS transitions produce a short GRB, leading to 

T^SGRB /„ -,„-2 , o ^, in-n f (/b) \ f (/sGRb) 



_ X 10'^ ^ 3 X 10-1) [^] [ ^^^^^ . (17) 

7?(km) V V 50 / V 0.1 / ^ ' 

'^NS-^QS,max \ / \ / 

This additional factor /sgrb can be related to internal processes, such as the capacity of the burning front to proceed 
to the crust (see 33). 

Note that the value of the efficiency /sgrb ~ 1 — 10 % has been obtained here assuming that all short GRBs are 
produced by NS — > QS transitions. However, we will show below ( TO.1|) that this scenario naturally leads to short 
GRB durations Tgo < 0.1 s, so that only very short GRBs are good candidates for counterparts of these NS QS 
transitions. Then, it reduces in principle even more the value of /sgrb- 

On the other hand, observations of afterglows of ultra-short GRBs, identifications of their host galaxy and measure- 
ments of their distance are extremely rare so that their intrinsic rate is unknown. The distribution o f the duration 
of BATSE bursts show that '--^ 8% of the short bursts have a duration below 100 ms (jHorvathll2002f ). This factor 
cannot however be directly applied to the estimate of the intrinsic rate given by Eq. [51 Indeed, this population of very 
short bursts may very well be a separate group with different proper ties: they are usua lly made of a single short and 
hard spike with possible substructure on a timescale of a few 10 fj,s (jCline et al.lll999() . The analysis ((F/Vmax) and 
log N — log S dia gram) of the distr i butio n of very short bursts observed by BATSE and Konus gives some evidence 
for a local origin (ICline et all 119991 12005D . This is however a bit contradictory with more recent results obtained by 
Swift: 10 very short bursts with a duration of less 100 ms have been detected by Swift until June 201^]. All of these 
show indeed a single short duration hard spike in the BAT. In many cases, the afterglow has not been identified or 
is very weak at the limit of the detection (GRB 050925, GRB 051105A, GRB 070209, GRB 070810B, GRB 070923, 
GRB 090417A, GRB 100628A, GRB 120305A). There are however two very short bursts where the afterglow has been 
well detected and localized and where there is a good candidate for the host galaxy. In both cases, the host candidate 



^ Source: Swift GRB table at http: \\swif t .gsf c .nasa.gov\docs\swif t\archlve\grb_table .html 
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is an early type galaxy and the afterglow shows a large offset of a few 10 kpc: GRB 050509B seems to be associated 
to an eUiptical galaxy at z = 0.225 with an offset of 35-55 kpc (jGehrels et al. 1 120051: lBergeHl2011[): G RB 090515 is 
probably associated to an early type galaxy at z = 0.403 with a large offset of 75 kpc (lBergeHl20liT ). These two 
examples indicate a population at cosmological distance which is uncorrclated to star formation, in agreement with 
the discussion above. The sample is however much too small to allow for a determination of the intrinsic rate of very 
short GRBs and an estimate of /sgrb- 

5. CENTRAL ENGINE AND CRUST MASSES 

In this section, we model the basis of the central engine mechanism of the internal burning front that will induce a 
resultant outflow w i th several relati v istically mov i ng emitting regions. This di s cussion is based on a series of work s 
(|Cheng et all fl9M . iDai et all ITool lOuved et all [200l iPaczvnski et~all [20051 iXu et al.l[2009t iFischer et afllMol ). 
Nuclear processes involving quark deconfinement in the hadronic phase may happen if the ^ 1 ~ 100 GeV/c^ 
DM particle candidate self-annihilates liberating ^ MeV — GeV photons and other light particle pair products. 
The corresponding lump of strange quark matter (uds matter) evolves towards the formation of a fireball that may 
cause stress and tension on the base of the inner crust of a NS. Since there are only a few preliminar simulations 
(jAbdikamalov et al]|2009t iHerzog et al.ll201ll ) developed so far on the full process of energy transport from the inner 
deconfined regions to outer regions through a burning front there is not much information where this front may stop 
or whether it fully proceeds to the outer crust. In these simulations the NS burning condition is dynamically analyzed 
showing the possibility that a hybrid star may form if the conversion front is not able to proceed as burning. However, 
a full detailed treatment has not hitherto been performed. 

In previous work (jPerez-Garcia fc Silkll201ll ). it was found that the seeding is most likely a non central process, and 
th erefore, geometrical ly, the burning front progression may not be a radially symmetric dynamical process (see Fig. 2 
in iPerez-Garcial ()2012D ). The possibility of a beamed ejection of the outer crust aris es then from th e aniso tropy in the 
progression of the burning front. This has been somewhat explored in the work of iLugones et al.l ()2002[ ) where they 
discuss the possibility of preferred ejection of a fireball through the polar caps. In the present work, we discuss the 
constraints on the ejected mass as a result of the NS— >'QS conversion. Similar to what may happen in the heavy ion 
collision events in large colliders, a fireball may be formed and grow rapidly. Then, a pure nucleonic EoS of matter, 
would no longer be valid as nucleon quark content may be deconfined. Since the newer EoS is softened, later evolution 
may lead to the fact that the original hydrostatic structure is not longer energetically possible and re-adjustement of 
the object to favor lower radii and masses and then to build up tension in the crust tending to eject it and break it up. 

The crust mainly corresponds to matter where the density is below the NSD. Therefore for lower densities a myriad of 
nuclei with mass number A populate this phase. Since the particular distribution of nuclei depend on the competition 
of short-range hadronic interaction and long range Coulomb interaction a set of irregular shapes different from the 
spherical one can arise, forming what is known as the pasta phase. As a result, low density matter in the crust of NSs 
is mostly neutron-rich due to the deleptonization caused by neutrino escape in the first stages of cooling. Its isospin 
content is closer to that of neutron matter than to regular nuclei where proton and neutron content are balanced 
Z ^ N. From the ^^Fe iron content usually assumed in the lower density region in the crust, there is a series of 
neutron-rich nuclei going through heavier ^^Ni, ^^Ge, ^^°Sr and after this the ND transition takes place, signaling the 
inner crust of the NS. In the inner crust, surrounding th ese non-uniform pasta structures, a neutron ga s is filling the 
system, as has been directly simulated in previous work (jHorowitz et al.ir2004 120051 lPerez-Garciall2006l) . 

Fragmented emission of ejecta is therefore possible since the stress in the base of the crust from a burning front may 
liberate the most abundant structures and the rest of lower A nuclei. The gradient of composition (heavy — > light) 
should also produce some variability during the crust ejection, all regions of the crust being not necessarily ejected 
with the same Lorentz factor. An additional source of uncertainty is related to the fraction of the initial energy release 
associated to the transition in the core which will be injected into gamma-ray photons, and the associated photo- 
disintegration of heavy nuclei. An efficient photo-disintegration would lead to a modified chemical composition biased 
towards light elements, and would enrich the medium by free neutrons that can be an additional source of energy 
(neutron decay) and internal dissipation in the ejecta (see §6.ip . In order to size the importance of this ejection, we 
plot in Fig. [3] the crust mass as a function of the stellar radius. These d ata are obtained in tegrating the TOV 
(jOppenheimer &: Volkofj[T939[ ) equations for eight representative different EoS (jPatta et al.lll995l ) up to the core. We 
compute the mass of the crust by integrating above a critical density, either the neutron drip density pnd ~ 4 10^^ 
g/cm'^ to estimate the mass of the outer (solid) crust only, or the nuclear saturation density pnsd ~ 2 x 10^^ g/cm'^ 
(where crust-core boundary sets in) to estimate the mass of the entire crust. We find that the condition M^j < lO"** M© 
needed for the relativistic motion (S|3|) can be fulfilled if only the outer crust is ejected {Mc ^ lO^^M© — 10^^'^ M©), 
whereas it becomes much more difficuh if the whole crust is expelled (Mc ~ lO-^-^M© - lO-^-^M©). The details 
of the physics of the propagation of the burning front and the associated energy deposition are required to estimate 
precisely the fraction of the crust which is ejected. If this fraction varies from a NS^QS transition to another, this 
would naturally lead to /sgrb < 1 as discussed in the previous section. 

6. GRB DURATION, LIGHT CURVE, SPECTRUM 

In this section we discuss some expected properties for SGRBs produced in the scenario proposed in this work. We 
focus on the identification of distinct features related to the specific central engine that could help in distinguishing 
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EoS considered in lDatta et al.l I I1995I) as a function of stellar radius 



this mechanism from others proposed in the hterature. 

6.1. Duration and light curve 

Let us assume that the ejected outer crust has a mass M^j = x 10~^ Mq, an initial energy Ecj ~ /cj,-3 x 

3.5 X 10^° erg and a width A = cAi, where At = At-Q x 10~^s is the hght crossing time (A = At-Q x 300m). The 
maximum final Lorentz factor that can be reached in the ejection is 

Due to the large physical uncertainties regarding the energy deposition in the crust and the following ejection, we do 
not attempt here a detailed calculation of the relativistic ejection (see e.g. iPan fc Sarill2006l for a self-similar solution 
of the propagation of a strong shock wave within the NS and the following shock breakout). For a thermal acceleration, 
the saturation to this value of the Lorentz factor will occur at radius 

i?,at ~ Pi? ~ 2 x 10^ ^47,-5/cj,-3 cm , (19) 

assuming an initial radius i?cj = i? — A ~ i?NS = 12 km for the ejection. The ejccta will become transparent for its 
own radiation at the photospheric radius 

i?p,c.y^c.2xl0i3M^/!5Cm, (20) 

where we take the Thomson opacity k ~ 0.2 cm^/g. Here we assume = 0.5 free electron per nucleon in the expanding 
gas, whereas the dynamical chemical composition discussed in the previous section may lead to lower values. This 
does not affect too much the discussion as the dependence goes moderately as i?ph oc k^-^. 

This expression of Rph is valid for i?ph ^ Ris — 2r^A. The radius Ris is defined below and is of the order of a few 
10^ cm so that the condition is always true for the typical parameters considered here. The estimates of i?sat and i?ph 
are based on standard fireball theory for GRBs (see e.g. iPiranI [2"004[ ) and are only rough estimates for the scenario 
considered here as t he physics of the e jection of the crust is much more complex. We should bear in mind that the 
ejecta composition (iKotera et al.ll2012D a nd fragmentation may affect both the dynamics (i?sat) and the interaction 
with radiation (Rph) (jOuved et al.l 120051 ). Therefore, if the energy is not deposited in an homogeneous way in the 
expelled crust, the final Lorentz factor in the ejecta may not be uniform. In addition, fragmentation of the crust during 
its ejection will also lead to some variability in the ejecta. If variations of P arc present on length scales civar < A, 
this will induce collisions (internal shocks) that will dissipate energy at a typical radius Ris given by 

< 2P2A ~ 2 X 10" Af T2_g/2. _3At_6 cm . (21) 

Clearly, this possible internal dissipation will always occur well below the photosphere, and, depending on the value 
of the Lorentz factor, even before the acceleration is complete. Most of the dissipated energy should contribute again 
to the acceleration and the internal shock phase will only tend to smooth out the initial internal variability, without 
contributing to the emission. 

The ejecta will expand freely initially, but will eventually be decelerated by the external medium. However, due to 
the expected delay between the formation of the NS and the NS — ?• QS transition, this external medium can correspond 
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to the periphery of the host galaxy or even the surrounding intergalactic medium, i.e. have a low density. We assume 
a uniform medium with a number density n ~ n_3 x 10^'^ cm"'^. Then the deceleration will start at radius 

- l»"M3?=/i'l":3'=om^ ,22, 

TUp is the proton mass. From these different estimates, we find for the proposed scenario that 

-Rsat ^ ^is ^ -Rph ^ -Rdcc ■ (23) 

This would naturally lead to two episodes of emission, a single spike emitted at the photosphere followed by an afterglow 
starting at late time due to the high value of i?dec- The duration of the prompt spike should be fixed by the intrinsic 
curvature of the emitting region and its lateral extension, which gives 



min Ml_J.l,; ( -i ) x 0.8 m'% s . (24) 



Except if the ejection is highly beamed, the minimum is usually given by the first term. Nevertheless, this estimate 
clearly points out towards short (< 1 s) and even probably very short (< 100 ms, see Figure |4]) GRBs without any 
strong variability (one main single spike). Due to the low external density, the afterglow should rise slowly and reach 

a maximum around tdcc = ^dcc/2r^c ~ 2 x 10"* A/^j^i^g/^^^gnlg^'^ s, i.e. a few 10 ms to a few hours after the GRB, 
depending on the external density. In addition, the combination of a low external density and a moderate energy 
will naturally lead to a weak afterglow. Note that if free neutrons are present in the ejecta, due to efficient photo- 
disintegration (see fJS]), these neutrons will decay at a typical radius R ~ Fcr^ ~ 5.4 x lO^"' A/^y^„g/cj._3 cm < i?dcc, 
fixed by the mean lifetime Tp ~ 900 s, which could lead to an early additional signature at iobs — 20 A/cj._5/j,^^_3 s (^^^ 
e.g. lBeloborodov|[200l . 

6.2. Spectrum 

If most of the gamma-rays are produced at the pho tosphere, one should expect a thermal (quasi-Planckian) spec trum, 
possibly modified at high energy by comptonization (jRees fc Meszaroj|2005l : |Pe'er et al.ll2006HBeloborodovll2011l ). The 
peak energy of the observed spectrum will then be located at Ep ~ S.OfcBTph, where Tph is the temperature of the 
photosphere. It can be computed assuming an adiabatic radial expansion from the ejection to the photosphere: 

- ilA) "i^f"-'' X»V ^ (25) 

The efficiency of the photospheric emission can also be deduced from the adiabatic evolution and equals ~ 
(^ph/-Rsat) This leads to a better estimate of the isotropic equivalent energy radiated in gamma-rays, 

^ 2 X 104« erg/s . (26) 

that stands on the lower values for GRBs. Of course, there is still a large uncertainty due to the unknown factors 
/b, /cj and Mej. However, these two estimates confirm the capacity to produce bright and hard spikes of gamma-rays 
in the considered scenario, as illustrated in Figure U] where all the constraints on the prompt gamma-ray emission 
expected in the NS — >■ QS conversion scenario are summarized. 

Another possible signature of the scenario would be the presence of spectral features an d lines associat ed to the 
specific chemical composition (heavy elements) of the ejected material. Following the work bv lMeszaros fc Reea (,1991 
where F ~ 10 — 100 were considered, the spectrum could be influenced below the MeV range. 

6.3. Gravitational waves and other non-photonic signatures 

In addition to the expected short duration burst of 7-rays and X-rays and to the afterglow, it is likely that a multi- 
messenger approach must be followed to spot this progenitor scenario for SGRBs. The emission of gravitational waves 
by the merger of NS-I-NS or NS-I-BH binary is the most promisin g source for detectors such as Virgo and LIG O and 
their advanced versions (jAcernese et all 120091: iHarrv et al.ll2010[) . or the Einstein telescope (|Hild et all [20100 . The 
predicted ern ission during the three main stages of the merger has been extensively studied in the literature (see 
lHughejl2009l for a review). In the inspiral stage, the signal is very well known up to the last stable orbit. Then, the 
emission during the merger is more uncertain a nd must be modeled using supercom puter simulations which provide 
information about the grav i tational waveforrn (iDuez et al.l 120061: iBaiotti et al.l 120101) fsee the rec e nt simulat i ons b y 
Shibata fc Taniguchil (|2006|): iPrice fc Rosswod (120061): iBaiotti et all (I2008D: [Anderson et all (I2OO8I ): iLiu et al.l (|2008l) 



Giacomazzo et al.l ( 2009[ ): iRezzoUa et al. I (120101 ): IKiuchi et al.l (120101) : IG iacomazzo et al.l (|2011[ ) which compute the 
evolution of the merger up to black hole formation, including or not magnetic fields, realistic EOS, etc.). The final 
ringdown stage is also very well known as the inspiral. Therefore, the detection of gravitational waves in association 
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Fig. 4. — In the parameter space mass of the ejected outer crust A/ej vs efficiency of the energy injection in the crust /cj, the following 
constraints are plotted; (1) lines of constant Lorentz factor are plotted in blue for F = 1 (non relativistic limit), 10, 100 and 1000. The 
limit F ~ 15 obtained from the compactness argument (see JSjis plotted in thick line and the forbidden region is shaded; (2) the limit where 
the radius of the internal dissipation ijjg is equal to the radius of the photosphere is plotted in red (most of the parameter space is well 
below this line, i.e. R^g <S ^?ph); the limit where the observed duration of the prompt emission of photospheric origin is equal to 100 ms (1 
s) is plotted in magenta solid (dashed) line; the limit where the observed peak energy of the photospheric emission is 100 keV is plotted in 
green; the limit where the isotropic equivalent gamma-ray energy emitted by the photosphere is equal to lO^^erg.s"^ is plotted in black. 
The effect of the redshift of the source on the duration and the spectrum are not included. Other parameters are At-g = 1 (i.e. the width 
of the outer crust is 300 m) and = 50 (i.e. the ejecta is beamed within ~ 10°). As can be observed, a large fraction of the parameter 
space (top-left region) corresponds to ultra-relativistic outflows (F > 100) producing a very short (< 100 ms) but bright {E-y iso > 10^^ erg) 
spike of hard {Ep > 100 keV) photons, i.e. a very short GRB. 

with a SGRB would undoubtedly prove if the progenitor is a merger or not. With a horizon of ^ 200 Mp c for NS-I-NS 
merg ers and ^ 420 Mpc for NS-I-BH mergers, there may be merger detections with advanced Vir go, LI GO ()Harrv et al.l 
12010( 1 (0.4 to 400 mcrgcrs/yr for NS+NS, a slightly smaller rate for NS+BH, lAbadie et al]l2010f l. However, the 
simultaneous detection with a SGRB if the merger scenario is correct is more uncertain, both for theoretical (beaming 
of the gamma-ray emission) and instrumental (localization of SGRBs) reasons. 

The gravitational wave signature of the scenario studied here is not as well known as for compact binary mergers. 
It is expected that a change in the moment of inertia is caused by deformations in the NS— t'QS transition. Some 
preliminary estimates of the transient gravitational wave signal from an explosive quark-hadron phase transition have 
been done (Staff et al. 2012). It could in principle be detected out to 20 Mpc with advanced Virgo/LIGO, which makes 
unfortunately the probability of detection with an associated SGRB quite low, even if the intrinsic rate of very short 
GRBs is quite uncertain (see SI)- 

The strategy for the detection of gravitational waves associated to a GRB is to combine searches, already having 
identified typical patterns in temporal and direction al coincidence with SGRBs that had sufficient gravitational- wave 
data available, although it is a very challenging task (lAbadie et al.|[2010l ). 

An additional signal that may be used to discriminate between progenitor models for SGRBs is the neutrino emis- 
sion. The specific signal expected from NS+NS mergers has been studied using supercomputer simulations (see e.g. 
iDessart et ani2009[ ). In this work and within the present scenario we do not attempt to describe such a neutrino emis- 
sion since it crucially depends on the central engine details. However, it should be expected a neutrino flux originating 
from hadronic reactions happening in the central engine (conversion in the nucleon burning front of quark deconfine- 
ment and strangeness production) or in the ejecta (photo-meson interactions between shock-accelerated protons/nuclei 
and gamma-ray photons). 

In addition, if produced neutrons do not interact they will decay being accompanied b y anti-neutrinos . Em ission 
of prompt muon neutrino fluxes have been performed in some GRB general scenarios (jBaerwald et al.l 120111 ) and 
seem to be testable as recent preliminary estimated sensitivities for Ice Cube 86-strings quote values of E'^^{E) ^ 

^ Source:http://www.ligo.org/science/GW-Burst.php 
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5 X l O^^GeV cm^^ sr~^. Additionally, experimen ts as AMANDA-II an d IceCube have the capability (iHughev et alj 
I2007D of detecting anisotropies from the emission (jAartsen et al.l I2012D of neutrinos from gamma ray induced air 
showers. Therefore, as motivated, the discovery of high energy neutrinos in correlation with a SGRB and GW signal 
would help disentangling the current scenario. 

7. CONCLUSIONS 

In this work we discuss the possibility that a NS— )-QS transition may be a central engine model for short (< 1 s) 
or more probably very short (< 100 ms) GRBs. This is an alternative to the popular NS+NS or NS+BH merger 
scenario for SGRBs. Note that this scenario is known to be compatible with observations of SGRBs but has not been 
proved yet, and that in addition, there are a few observations suggesting that very short bursts could be a separate 
population. We suggest that these very short bursts could be produced by the ejection and acceleration to relativistic 
speed of the (outer) crust of the NS during the conversion to a QS. Wc find that 

• The isotropic equivalent gamma-ray energy ^ 10"*^ — 10^^ can be accounted for, assuming that 0.1 — 1% of the 
gravitational energy released by the transition is injected in the outer crust with, however, a large uncertainty 
due to the unknown beaming factor of the ejection; 

• High Lorentz factors, necessary for the emission of 7-rays on short time-scales (compactness problem), can be 
reached, as long as the mass of the expelled crust is less than ^ IO^^A/q; this is attainable for the outer crust 
in NS models. 

• The rate of very short GRBs can probably be reproduced, assuming that only a fraction of transitions lead to 
a GRB and that there is on average a large delay between the formation of the NS and the conversion into a 
QS. Such a delay is expected due to the fact that the transition is triggered by self-annihilation of DM, which 
needs first to be accreted in the core of the NS. Such delays would then naturally at least partially suppress the 
correlation between star formation and very short GRBs, which should be observed in any type of galaxy with a 
large offset (at the present time, there arc only two host galaxies of very short GRBs which have been possibly 
identified, both being early-type, and the offset is a few 10 kpc in both cases); 

• The prompt gamma-ray emission should be mainly produced at the photosphere, without a strong variability 
which should be washed out by internal dissipation at much smaller radii. 

• For a large fraction of the parameter space, a hard {Ep > 100 keV) and short (duration < 100 ms) spike of 
gamma-rays is expected, in general agreement with observations of very short GRBs; 

• The afterglow should rise slowly and be rather weak, due to a low density of the external medium. A possible 
additional signature can be expected at early time if the ejected material initially contains free neutrons. 

• Fragmentation of the ejecta arises naturally in this model since non-uniform nuclei arranged in lattice or even 
pasta phases are present in the low density matter. 

• Spectral features of heavy nuclei relativistically accelerated are expected below ~ 1 MeV. 

To summarize, possible signatures compared to the binary (NS, BH) merger scenario are the shortness of the 
prompt gamma-ray emission, with possibly a thermal spectrum and spectral features due to the heavy composition, 
the associated GW emission, and possibly the properties of the host galaxies and the distribution of the afterglow 
position in the host, this latter signature being however difficult to characterize due to the uncertainties on the typical 
delay between the NS formation and the transition to a QS. Clearly, a multi-wavelength/multi-messenger approach is 
needed to reach a firm conclusion. 

We thank the MICINN (Spain) MULTIDARK, FIS-2009-07238, FIS2011-14759-E and FIS2012-30926 projects and 
the ESF-funded COMPSTAR project for partial financial support. We acknowledge useful discussions with I. Bombaci, 
J. Horvath, K. Kotera, R. Ouyed and C. Providencia. M. A. P. G. would like to thank the lAP for its kind hospitality. 

REFERENCES 



Aalseth, C. E. ot al., CoGeNT Collaboration 2011, Phys. Rev. 

Lett. 106, 13, 1301 
Aartscn, M. G. et al., IceCube Collaboration 2012, arXiv: 

1210.7992 

Abadie, J. et al. 2010 , ApJ, 715, 1453 

Abadie, J., Abbott, B. P., Abbott, R. et al. 2010, Classical and 

Quantum Gravity, 27, 173001 
Abdikamalov, E. B., Dimmelmeier, H., RezzoUa, L., and Miller, J. 

C. 2009, MNRAS, 392, 52 
Acernese, F. et al. (Virgo Collaboration) (2009), Virgo internal 

note VIR-0027A-09, https://tds.ego-gw.it/itf/tds/ 



Ackermann, M., Asano, K., Atwood, W. B., et al. 2010, ApJ716, 
1178 

Alcock, C, Farhi, E., & Olinto, A. 1986, ApJ, 310, 261 
Anderson, M. et al., Phys. Rev. Lett. 100, 191101 
Arzoumanian, Z. et al. 2002, ApJ, 568, 289 
Avrorin, A. et al. 2009, Proceedings of the 31st ICRC, Lodz 
Baer, H. et al. 2004, JCAP 0408, 005 
Baerwald, P., et al. 2011, Phys. Rev. D 83, 067303 
Baiotti, L., Giacomazzo, B., & RezzoUa, L. 2008, Phys. Rev. D, 
78, 084033 



13 



Baiotti, L., Damour, T., Giacomazzo, B., Nagar, A., & RczzoUa, 

L. 2010, Physical Review Letters, 105, 261101 
Baring, M. G., & Harding, A. K. 1997, ApJ,491, 663 
Beane, S. R. et al. 2011, Phys. Rev. Lett. 106, 162001 
Bchnke, E. et al. 2012, arXiv:1204.3094 .^1 
Beloborodov, A. M. 2003, ApJ, 585, L19 
Beloborodov, A. M. 2011, ApJ, 737, 68 
Berezhiani, Z. et al. 2003, ApJ, 586, 1250 
Berger, E. 2007, ApJ670, 1254 
Berger, E. 2011, New Astronomy Reviews, 55, 1 
DAMA Collaboration, Bernabei R. et al. 2010, Eur. Phys. J. C 

67, 39 

Bertone, G., cd.. Particle Dark Matter; Observations, Models and 
Searches. Cambridge University Press, 2010, ISBN 
978-0-521-76368-4. 

Bombaci, I., Parenti, I., Vidaiia, I. 2004, ApJ614 

Bombaci, I. and Popov, S. B. 2004, A&A, 424, 627 

Bombaci, I. and Datta, B. 2000, ApJL, 530, 69 

Bringmann, T. et al. 2012, arXiv: 1203. 1312 

Bromberg, O., Nakar, E., Piran, T., & Sari, R. 2012, ApJ, 749, 
110 

Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 
1111 

Cheng, K. S., Dai, Z. G., 1996, Phys. Rev. Lett. 77, 1210 
Chodos, A. et al. 1974, Phys Rev. D 10 2599 
Cline, D. B., Matthey, C, & Otwinowski, S. 1999, ApJ527, 827 
Cline, D. B., Czerny, B., Matthey, C, Janiuk, A., & Otwinowski, 

S. 2005, ApJ, 633, L73 
Dai, Z. G., T. Lu, T. 1998 , Physical Review Letters 81, 4301 
Daigne, F., & Mochkovitch, R. 1998, MNRAS, 296, 275 
Datta, B. et al. 1995, Astrophys. Astr. 16, 375 
Dahlen, T., Strolger, L.-G., Riess, A. G., et al. 2004, ApJ613, 189 
Demorest, P. et al. 2010, Nature 467, 1081 

Dessart, L., Ott, C. D., Burrows, A., Rosswog, S., & Livne, E. 

2009, ApJ, 690, 1681 
Duez, M. D., et al. 2006, Phys. Rev. Lett. 96, 031101 
Fang K., Kotera K., Olinto A. V. 2012, arXiv:1201. 5197^^1, 

Kotera, K., Perez-Garcia, M. A., Silk, J. arXiv: 1303.1186 
Fischer, T. et al. 2010, Class. Quantum Grav. 27, 114102 
Frey, R. E. et al. 2007, AIP Conf. Proc. 928, 11 
Gehrels, N. , Ramirez-Ruiz, E. and Fox, D.B. 2009, ARAA, 47, 

567 

Gehrels, N. et al. 2005, Nature 437, 851 

Ghirlanda, G., Ghisellini, G., & Celotti, A. 2004, A&A, 422, L55 
Giacomazzo, B., Rezzolla, L., & Baiotti, L. 2009, MNRAS, 399, 
L164 

Giacomazzo, B., RczzoUa, L., & Baiotti, L. 2011, Phys. Rev. D, 
83, 044014 

Glendenning, N.K. 2000, Compact stars, Ed. Springer- Verlag, 
New York. 

Glendenning, N.K., Moszkowski, S. A. 1991, Phys. Rev. Lett. 67, 
2414 

Goldman, I. and Nussinov, S. 1989, Phys. Rev. D 40, 3221 
Gould, A. 1987, ApJ, 321, 571 
Guetta, D. and Piran, T. 2006, A&A, 453, 823 
Guiriec, S., Briggs, M. S., Connaugthon, V., et al. 2010, ApJ, 
725, 225 

Haensel, P., Potekhin, A.Y., Yakovlev, D.G. 2007, Neutron Stars 
I: Equation of State and Structure, Ed. Springer- Verlag, New 
York. 

Harry, G. M. and the LIGO Scientific Collaboration 2010, Class. 

Quant. Grav., 27, 084006 
Hascoet, R., Daigne, F., Mochkovitch, R., & Vennin, V. 2012, 

MNRAS, 421, 525 
Herzog, M. and Roepke, K. 2011, 'arXiv:1109.0539Vl 
Hild, S., Chelkowski, S., Freise, A., Franc, J., Morgado, N., 

Flaminio, R., & DeSalvo, R. 2010, Class. Quant. Grav., 27, 

015003 

Horvath, J. E. 2002, A&A, 392, 791 

Horowitz, C. et al. 2005, Phys. Rev. C 72, 035801 

Horowitz, C. et al. Phys. Rev. C 69, 045804 

Horvath J. E. , Benvenuto, O. G., Vucetich, H., Physical Review 
D 45, 10, 3865 

Hughey, B. for the IceCube Collaboration 2007, Journal of 

Physics: Conference Series 60, 340 
Hughes, S. A. 2009, ARA&A, 47, 107 



loka, K. and Taniguchi, K. 2000,ApJJ. 537, 327 
Jaikumar, P., Rapp, R., and Zahed, L 2002, PRC, 65, 055205 
Kaplan, D. E. , Luty, M. A. and Zure,k K. M. 2009, Phys. Rev. 
D. 79, 115016 

Kiuchi, K., Sekiguchi, Y., Shibata, M., & Taniguchi, K. 2010, 

Physical Review Letters, 104, 141101 
Kobayashi, S., Piran, T., & Sari, R. 1997, A pJ, 490, 9 2 
Kouvaris, C. and Tinyakov, P. 2010, arXiv: 1004.05861 /1 
Kouvaris, C. 2008, Phys. Rev. D 77, 023006 
KouveUotou, C. et al. 1993, ApJ, 413, LlOl 
Kuhlen, M. 2010, Advances in Astronomy, 2010, 162083 
Landau, L. D. 1932, Phys. Z. Sowjetunion, 1, 285. 
Lattimer, J. M. and Prakash, M. 2005, Phys. Rev. Lett. 94, 

111101 

Lionel, B. and Tinyakov, P. 2012, Physical Review Letters, 109, 6, 
061301 

Lithwick, Y., & Sari, R. 2001, ApJ, 555, 540 

Liu, Y. T., Shapiro, S. L., Etienne, Z. B., & Taniguchi, K. 2008, 

Phys. Rev. D, 78, 024012 
Lopes, I. and Silk, J. 2010, 330, 462 

Lorimer, D. R. 2004, in Young Neutron Stars and Their 
Environments, lAU Symposium, Vol. 218, Camilo, F. and 
Gaensler B. M., eds. 

Lugones, G. et al. 2002, ApJL, 581, 101 

Madsen, J. 2001, Phys. Rev. Lett. 87, 172003 

Meszaros, P. and Rees, M. J. 1998, ApJL, 502, 105 

Mctzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, N., 
& Quataert, E. 2011, MNRAS, 413, 2031 

Nakar, E., Gal- Yam, A., & Fox, D. B. 2006, ApJ, 650, 281 

Nakar, E. 2007, Phys. Rep. 442, 166 

Navarro, J. et al. 2004 , MNRAS, 349, 1039 

Oppenheimer, J. R. and Volkoff, G. M. 1939, Phys Rev. 55, 374 

Ouyed, R., Dey, J., Dey J. 2002, A&A, 390, 39 

Ouyed, R., Rapp, R., Vogt, C. 2005, ApJ632, 1001 

Ouyed R. et al. 2005, ApJ, 626, 389 

Paczynski, B., Haensel, P. 2005, MNRAS, 362, 4 

Page, D. et al 2009, ApJ707, 1131 

Pan, M. and Sari, R., ApJ643, 416 

Pe'er, A., Meszaros, P., & Rees, M. J. 2006, ApJ, 642, 995 
Perez-Garcia, M. A. 2006, J. of Mathematical Chemistry 40, 63 
Perez-Garcia, M. A., Silk, J., & Stone, J. R. 2010, Physical 

Review Letters, 105, 141101 
Perez-Garcia, M. A. 2012, in the Proceedings of the Rencontres 

de Moriond Cosmology 2012, La Thuile, Aosta valley, Italy, 

arXiv: 120 5.25811 / 1 [astro-ph.HE] 
Perez-Garcia, M. A., Silk, J. 2012, Phys. Lett. B, 711 6 
Piran, T. 2004, Reviews of Modern Physics, 76, 1143 
Price, D. J. and Rosswog, S., Science 312, 719 
Rees, M. J. 1966, Nature, 211, 468 
Rees, M. J., & Meszaros, P. 1994, ApJL, 430, 93 
Rees, M. J., & Meszaros, P. 2005, ApJ, 628, 847 
Rezzolla, L., Baiotti, L., Giacomazzo, B., Link, D., & Font, J. A. 

2010, Classical and Quantum Gravity, 27, 114105 
Rhoads, J. E. 1997, ApJ, 487, LI 

Shibata, M., & Taniguchi, K. 2006, Phys. Rev. D, 73, 064027 
Soderberg, A. M., Berger, E., Kasliwal, M., et al. 2006, ApJ, 650, 
261 

Staff J. E., Jaikumar P., Chan V., and Ouyed R. 2012, ApJ 751 24 
Stanek, K. Z. 2003, et al., ApJL, 591, 17 
Thompson, C. 1994, MNRAS, 270, 480 
Usov, V. V. 1992, Nature, 357, 472 

Vogt, C, Rapp, R., and Ouyed, R. 2004, Nuclear Physics A, 735, 
543. 

Villasenor, J. S., Lamb, D. Q., Ricker, G. R., et al. 2005, Nature, 
437, 855 

Weniger, C. 2012. llarXiv:1204.2797l 
Witten, E. 1984, Phys. Rev. D, 30, 272 
Woosley, S. E. 1993, ApJ, 405, 273 

Xu R., Liang, E., Science in China G: Physics and Astronomy 52, 
315 

Yakovlev, D.G. and Pethick, C.J. , Annual Review of Astronomy 

and Astrophysics, 42, 169 
Zhang, B., & Yan, H. 2011, ApJ, 726, 90 

Zhang, B., Zhang, B. B., VirgiU, F. J., et al. 2009, ApJ, 703, 1696 



